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A.  INTRODUCTION 


One  of  the  most  effective  ways  to  save  electrical  energy  is  to  turn  off 
or  dim  electric  lights  whenever  possible.  Properly  designed  daylighting  can 
supplement  or,  in  some  cases,  completely  replace  electric  lighting  during  the 
daylight  hours  while  maintaining  the  predetermined  necessary  level  of  illumination. 

Design  techniques  have  been  developed  to  calculate  the  horizontal  illuminance 
on  the  work  plane  in  a  room,  based  on  luminous  flux  being  furnished  by  vertical 
openings  (windows)  or  horizontal  openings  (skylights).  This  report  discusses 
the  procedures  taken  to  calculate  the  average  illuminance  upon  a  workplane  due  to 
sloped  openings  (i.e.  skylights  in  a  pitched  roof). 

The  sloped  skylight  problem  is  divided  into  two  sections: 

1)  determination  of  total  external  daylight  illuminance  upon 
a  sloped  surface. 

2)  determination  of  average  illuminance,  upon  an  interior 
workplane,  due  to  the  sloped  skylight.  . 

The  division  of  the  problem  in  this  way  mades  it  relatively  easy  to 
develop  a  simple  procedure  for  calculating  the  average  illuminance  on  a  work- 
plane  due  to  skylights  located  in  a  pitched  roof. 
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B.  DETERMINATION  OF  TOTAL  EXTERNAL  ILLUMINANCE  ON  A  SLOPEO  SURFACE 


To  calculate  the  external  illuminance  upon  a  sloped  surface  many  factors 
must  be  known,  such  as  where  the  sun  is  relative  to  the  surface's  position, 
cloud  distribution  of  the  sky  (i.e.  cloudy-clear),  time  of  the  year,  and  the 
angle  the  sloped  surface  makes  with  the  earth's  surface  (a  horizontal  plane). 

Several  solar  quantities  are  calculated  to  determine  the  sun's  position 
relative  to  the  surface's  position  on  earth.  The  cloud  distribution  of  the 
sky  is  picked  to  be  the  two  extreme  cases  of  completely  clear  sky  and 
completely  overcast  sky. 

Once  the  above  mentioned  variables  are  known  the  external  illuminance 
on  a  sloped  surface  can  be  determined  for  the  completely  overcast  sky  and  the 
completely  clear  sky  cases. 

1)  Determination  of  Solar  Quantities: 
a)  Solar  Time  (Reference  1) 

Solar  time  (ST)  is  defined  as  the  time  based  on  the  apparent  angular  motion 
of  the  sun  across  the  sky,  with  solar  noon  the  time  the  sun  crosses  the  meridian 
(longitude)  of  the  observer.  Solar  time  does  not  coincide  with  the  local  standard 
clock  time.  It  is  necessary  to  convert  standard  time  to  solar  time  by  using 
Equation  la. 

ST  =  Standard  Time  +  AM  (la) 

where  ST  =  solar  time 

aM  =  correction  to  solar  time  in  minutes. 

AM  is  given  in  Equation  lb. 

AM  =  4(LST-LL0C)  +  ET  (lb) 

where  L^j  =  Standard  Meridian  for  the  local  time  zone. 

See  Chart  2  of  Figure  1  for  a  listing  of  Standard 
Meridians  of  Time  Zones. 


*"LOC  =  longitude  (degrees  west)  of  site  location. 

ET  =  Equation  of  Time  (Reference  2).  The  function 


is  depicted  in  Chart  1  of  Figure  1  and  can 
be  calculated  by  Equation  lc. 

ET  =  9.87  sin2B  -  7.53  cosB  -  1.5  sinB  (lc) 

where  B  =  (n-81)/58  rad:  n  is  the  number  of  the  day  in  the  year 
and  a  partial  listing  of  n  is  given  by  Table  1. 

The  equation  of  time  may  be  read  directly  from  Chart  1  of  Figure  1  to  get 


a  crude  approximation. 

For  example,  on  March  21 

ET  is  about 

-8  minutes. 

Table  2  gives  the  approximate  ET's  for  the  two  equinoxes  and 

TABLE  1 

solstices. 

MONTH 

DAYS  IN  MONTH* 

n 

JAN 

31 

i 

FEB 

28 

31  +  i 

MARCH 

31 

59  +  i 

APR 

30 

90  +  i 

MAY 

31 

120  +  i 

JUNE 

30 

151  +  i 

JULY 

31 

181  +  i 

AUG 

31 

212  +  i 

SEPT 

30 

243  +  i 

OCT 

31 

273  +  i 

NOV 

30 

304  +  i 

DEC 

31 

334  +  i 

where  i 

is  the  day  of  tne  month  in 

question. 

★ 

This  table  is  for  the 

average  year  of  365  days. 

For  leap  year  n  should 

increased  by  1  for  the  months  after  FEB. 
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Day  number 


Month 


Standard  Meridian 

Zone  Position  (Longitude) 


+4 

Atlantic 

60°  W 

+5 

Eastern 

75°  W 

+6 

Central 

90°  W 

+7 

Mountain 

105°  W 

+8 

Pacific 

120°  W 

FIG.  1. 

CHART  I  THE  PERIODIC 

SHIFT  FOR  THE  SOLAR 

Chart  1 


Chart  2 

Standard  Meridians 
of  Time  Zones 


CHART  2  A  LIST  OF  THE  STANDARD  MERIDIANS  OF  TIME  ZONES 


TABLE  2 


APPROXIMATE  VALUES  OF  ET 


SEASON 

DATE 

ETj 

SPRING 

MARCH  21 

-8 

SUMMER 

JUNE  21 

-1 

FALL 

SEPTEMBER  21 

+8 

WINTER 

DECEMBER  21 

+  1 

b)  Solar  Altitude  Angle 

The  solar  altitude  angle  is  the  vertical  angle  of  the  sun  above  the 
horizon,  as  seen  in  Figure  2.  The  solar  altitude  (SA)  may  be  found  by  using 
Equation  2. 


SA  =  ARCSIN(sin6sin(L)  -  cos(L)cos6cos(ST  x  15°))  (2) 

where  L  =  latitude  of  the  site 

6  =  the  declination  angle,  which  is  defined  as  the  angular 
position  of  the  sun  at  solar  noon  with  respect  to  the 
plane  of  the  equator.  6  can  be  found  by  using  Equation  3. 


6  =  23.45°sin(360x^^-) 

where  n  =  the  number  of  the  day  in  the  year  (listed  in  Table  1) 


(3) 


c)  Solar  Azimuth  Angle 

The  solar  azimuth  angle  (SAZ)  is  defined  as  the  angular  displacement  from 
south  of  the  projection  of  the  line  from  sun  to  site  on  the  horizontal  plane,  as 
shown  in  Figure  2.  The  SAZ  may  be  calculated  using  Equation  4. 

SAZ  =  ARCSIN[(cos6sin(-l 5°xST)/cosSA)  (4) 
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TO  S>Oi  OoW'E 


Figure  2.  Angular  representations 
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2)  Illuminance  from  an  Overcast  Sky: 
a)  Overcast  Sky  Luminance 

The  International  Commission  on  Illumination  (CIE)  adopted  the  Moon  and 
Spencer  formula,  given  in  Equation  5,  to  represent  the  luminance  distribution 
of  the  completely  overcast  sky  (Reference  3). 

L  =  (Lz/3)  [1+2COS0]  (5) 

where  Lz  =  sky  zenith  luminance 

e  =  Angular  distance  from  zenith  to  sky  element  (see  Figure  2). 

The  sky  zenith  luminance  is  also  given  by  the  CIE  as  Equation  6. 

L  =  123  +  8600  sin(SA)  -dyde-la  (6) 

z  nr 


b)  Overcast  Sky  Illuminance 

The  illuminance  on  a  sloped  surface  under  completely  overcast  sky 
conditions  is  given  in  Equation  7  (see  Appendix  A  for  derivations  of  overcast 
sky  equations) 


where 


E0C  EH0C  C0Sa  +  EV0C  Sina 


(7) 


ehoc 

Evoc 


=  Horizontal  illuminance  for  overcast  skies 
=  Vertical  illuminance  for  overcast  skies 

=  angle  of  tilt  of  sloped  surface  with  respect  to  earth  surface 


The  completely  overcast  sky's  horizontal  illuminance  is  calculated  by  Equation  8. 
L. 


e  =  Ji  r— 

tH0C  3  l3 


tt/2 

l  (*-2<r)[ 


.2  .2 

sin  o^-sin  e 


-  -  j  (cos3e2  -  cos3©^] 


e  =  2  ’a,einc 


(8) 


[Angles  in  Radians] 

*The  sumnation  is  read  as:  0  is  summed  from  -a  to  j  by  0inc  increments 
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where  L2  =  zenith  luminance  for  overcast  skies 

e^.eg  =  upper  and  lower  bounds  from  the  zenith  (see  Figure  3)  to 
sky  element. 

<p'  =  Average  angular  distance  from  East  to  the  surface's  projected 

intersection  with  the  sky  sphere  (see  Figure  4). 

<p'  is  calculated  by  Equation  9  (see  Appendix  B) 

=  ARCSIN(CTN0AVGCTNa)  Radians  (9) 

where  *  Average  angular  distance  from  zenith  to  sky  element 


The  vertical  illuminance  component  of  the  completely  overcast  sky  is 
calculated  using  Equation  10. 


EV0C 


tt/2 

I 


cos$'[(6i-02)  +  ^  (sin301-sin3e2)  -  ^  (sin2e-|-sin2e2)]] 


inc 

[Angles  in  Radians] 


OO) 


A  listing  for  the  overcast  sky  illuminance  (EqC)  is  given  in  Table  3.  It 
is  listed  with  respect  to  time  of  day,  day  of  the  year,  latitude  of  site,  and 
slope  of  the  surface  in  question. 
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TABLE  3  OVERCAST  SKY  ILLUMINANCE  (Ft-c) 
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3) 

a) 


Illuminance  from  a  Clear  Sky: 

Solar  Illuminance  (Direct  Sunlight) 


The  direct  perpendicular  (normal  to  the  earth's  surface)  solar  illuminance 
is  listed  in  Figure  7-8  of  the  IES  Handbook,  6th  ed.  (Reference  4),  or  can  be 
calculated  by  Equation  11,  which  is  a  polynomial  "fit"  of  Figure  7-8. 

Epsi  =  -17.9  +  427 . 5(SA)  -  8.7(SA)2  +  0.086(SA)3  -  0.3306E-3(SA)4  (11) 

where  SA  =  Solar  altitude  in  degrees. 

The  horizontal  and  vertical  solar  illuminances  are  given  by  Equations  12a  and 
12b,  respectively. 


EHSI  =  EPSISin(SA) 
EVSI  =  Epsic°s(SA) 


(12a) 

(12b) 


The  total  solar  illuminance  on  a  sloped  surface  with  any  orientation  is 
given  by  Equation  13  (see  Appendix  C). 


ETSI  =  EHSI  COSa  +  EVSI  sina  cos(6'SA^)- 
where  EHSI  =  horizontal  solar  illuminance 
EVSI  =  vertica^  solar  illuminance 
a  =  angle  of  surface  tilt 

6  =  surface  orientation  with  respect  to  south.  A  counter¬ 

clockwise  direction  is  positive. 

SAZ  =  solar  azimuth  angle  (Equation  4) 


(13) 


b)  Clear  Sky  Luminance 

A  standardized  luminance  distribution  of  clear  skies  accepted  by  the  CIE 
is  based  on  the  Kittler  equation  (angles  in  radians)  (Reference  5). 
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(1  -e~  •  32/cose  x  /  91+10e"3y+.  45cos2y ) 

L  Lz  .247(.9l+10e'3'^Z"SA^  +  .45cosZY) 

where  L  *  sky  luminance 

Lz  *  zenith  luminance 

y  *  angular  distance  from  sun  to  sky  element 
e  =  angular  distance  from  zenith  to  sky  element 
SA  *  solar  altitude 


(14) 


For  the  clear  sky  condition,  the  zenith  luminance  of  Equation  15a  is  used. 
(Reference  6). 


Lz  =  507e2‘35SA(cd/mZ) 


(15a) 


The  angle  y  is  calculated  by  Equation  15b. 

Y  *  ARCC0S[sinSAcoseAVG-cosSAsin6AVGsinUAVG-SAZ)]  ( 

0l+02 

where  eAVG  =  average  angular  distance  from  zenith  to  sky  element  (  ^  )• 

(See  Figure  4). 

<f>  =  average  angular  distance  from  East  to  sky  elements  ground 

AVG 

projection  (- ^ — )(see  Figure  4). 


c)  Clear  Sky  Illuminance 

The  total  illuminance  for  a  completely  clear  sky  is  composed  of  a  sky 
luminance  component  and  a  solar  illuminance  component  as  given  by  Equation  16 

Ec  =  [£HCcosa  +  Evcsina]  +  ETSJ  O 

where  Eur  =  horizontal  illuminance  due  to  sky  luminance 
HC 

=  vertical  illuminance  due  to  sky  luminance 
Etsi  =  total  solar  illuminance  (Equation  13) 
a  =  tilt  of  the  surface 
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Ehc  is  given  in  Equation  17  while  Eur  is  calculated  in  Equation  18. 


VC 


2  2 

n/2  sin  0,-sin  e,  2n 
:HC  =  I  C(  2  ^  I  L(*,-*i)]  ~ 

e=0*einc  *=0-*inc 


2  2 

u/2  sin  6,-sin  8,  h-*'+b 

C( - f - L)  l 

2  -a,0inc  *=*'+&, *. 

*  ¥  ’Mnc 


(17) 


-vc 


v/2  00-Qi  sin29i -sin20o 

I  (-V1  + - 3 - ^ 

°-8lnc 


2it+B 

l  U  (cos*, -cos*,) 

L-.lo  a.  1  ^ 


4>=Tr+6 ,  - 


me 


n/2  0o-0i  sin201-sin29o 

i  <-H->  *  — l — -> 

9=0,0. 


inc 


ir+B 

l  L| (cos*,-cos*,) 

b  =  B, *. 

K  'Mnc 


tt/2  0o-0i  sin2e,-sin20o 

i  (-V-* — 1 - h 


0_  2  -a,0irtc 


TT  -  <f)  +B 

I  L | (cos*, -cos*,) 

k  =  A"4.0  A  '  ^ 


08) 


*=*'+B,* 


me 


where  L  =  clear  sky  luminance 

6^02  =  upper  and  lower  bounds  from  the  zenith  to  the  sky  element 
(see  Figure  4). 

*j.*2  =  upper  and  lower  bounds,  from  East  to  the  sky  element's  pro¬ 
jection  on  the  earth's  surface 

*'  =  average  angular  distance  from  East  to  the  ground  projection 

of  the  surface's  intersection  with  the  sky  sphere  (see  Figure  5) 
6  =  surface  orientation  with  respect  to  South.  A  counter¬ 

clockwise  direction  is  positive. 

The  clear  sky  illuminance  of  Equation  16  is  tabulated  in  Tables  4,  5,  and  6. 
Table  4  contains  Ec's  for  March  21  and  September  21  (Spring  &  Fall),  Table  5 
contains  Ec's  for  June  21  (Summer),  and  Table  6  contains  Ec's  for  December  21 
(Winter).  The  Ec's  are  tabulated  with  respect  to  the  time  of  day,  day  of  the 
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year,  latitude  of  the  site,  slope  of  the  surface,  and  surface  orientation. 

The  surface's  orientation,  which  is  the  direction  the  slope  faces,  is  given  in 
these  tables  as  the  main  compass  points  N,S,E  or  W.  Figure  6  gives  a  division 
of  the  compass  into  four  quadrants  of  N,S,E  or  W  to  be  used  for  rough  approxi¬ 
mations  of  surface  orientation. 
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TABLE  4  CLEAR  SKY  ILLUMINANCE  (Ft-C)  -  MAR  &  SEPT  21 


TABLE  5  CLEAR  SKY  ILLUMINANCE  (Ft-C)  -  JUNE  21 


4)  Examples: 

a)  Overcast  Sky  Example 

What  is  the  illuminance  on  a  15°  tilted  surface  in  Boston,  Ma.,  for  an 
overcast  sky  on  Dec.  21  at  2  pm  standard  time? 

First  find  solar  time:  Equation  1.  Need  site's  latitude  and  longitude,  the 
standard  meridian,  and  ET,  as  well  as  standard  time. 

Site  Lat.  =  42°' 

•  from  map 

Site  Long.  =  71 °J 

Std  Merd.  =  75°  (chart  2  of  Figure  1) 

ET  =  +1  min  (Table  2) 

ISL _±: 

Solar  time  =  14:00  +  [4(75-71)  +  1]  =  14:00  +  : 17 
Solar  time  =  14:17  or  2:17  pm 

Now  use  21  Dec  column  and  15°  slope  row  of  Table  3.  Interpolating  between  2  pm 
and  3  pm  to  obtain  a  value  for  2:17  pm  yields. 

Eqc  (2:17)  =  584.5  ft-c 

b)  Clear  Sky  Example  (1 ) 

What  is  the  exterior  horizontal  illuminance  for  Washington,  DC,  on  June  21 
at  10  am  Daylight  Savings  Time  under  clear  sky  conditions. 

Find  solar  time:  Use  Equation  1 
Standard  Time  =  10-1  -  9  am 
Site  Lat.  =  39°' 

■  from  map 

Site  Long.  =  71° 

Std  Merd.  =  75° 

ET  =  -1  min 

Solar  time  =  (9:00)  +  [4(75-77)  +  (-1)] 

:  =  9:00  +  [-9] 

Solar  time  =  8:51  am 
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With  double  interpolation  from  Table  5,  as  shown  below,  the  horizontal 
illuminance  (0°  slope)  is  found  to  be  7305  ft-c.  Note  that  surface  orientation 
is  not  a  factor  in  the  horizontal  illuminance  case. 

i  Solar  Time 


Lat 

8:00 

8:51 

10:00 

38 

5878  + 

9277  + 

39 

5880  - 

7305 

-  9230 

42 

5886+ 

1 

9090+ 

Clear  Sky  Example  (2) 

Repeat  Example  (1)  but  with  a  15°  sloped  surface  with  a  northern  exposure. 
Use  Table  5  but  with  a  15°  slope  and  a  north  orientation.  From  before,  the 
solar  time  is  8:51  am  and  latitude  is  39°.  Using  double  interpolation  we  find 
an  exterior  illuminance  of  6841  ft-c. 


North  at  15° 


Solar  Time 


Lat 

8:00 

8:51 

10:00 

38 

5635  + 

8540+ 

39 

5647 

6841 

-  8457 

42 

5684  + 

8208  + 

Note  that  for  this  case  that  the  orientation  is  very  important,  it  determines 
which  column  to  use  in  Table  5. 
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C.  INTERIOR  ILLUMINANCE  FROM  SLOPED  SKYLIGHTS 


This  section  is  composed  of  two  parts.  First,  a  calculation  procedure, 
similar  to  that  in  Navy  Tech  Data  Sheet  80-09,  but  for  pitched,  rather  than 
horizontal,  skylights,  will  be  presented.  Second,  a  simplifying  straight  line 
equation  for  calculating  room  coefficients  of  utilization  will  be  developed. 

The  geometry  to  be  considered  is  shown  in  Figure  6.  It  consists  of  a  room 
of  length  £,  width  w  and  wall  height  h,  with  a  ceiling  sloped  at  an  angle  <j>. 
Skylights  of  length  a  and  width  b  are  imbedded  in  the  roof,  n  per  side.  It  is 
assumed  that  these  skylights  are  perfectly  diffusing  and  have  a  transmittance 
t,  a  reflectance  from  their  interior  side  ps  and  a  dirt  depreciation  factor 
SDD.  A  procedure  to  determine  the  average  maintained  illuminance  on  a  work 
plane  2.5  feet  above  the  floor  is  desired,  assuming  the  illuminance  on  the 
exterior  of  the  skylights  is  known. 

The  procedure  will  be  presented  through  three  work  sheets.  Worksheet  1 
deals  with  the  skylights  and  asks  the  user  to  enter  the  data  described  in  the 
previous  paragraph.  Worksheet  2  concerns  the  room.  The  first  six  items  are 
data  entries.  Items  7-10  involve  finding  the  effective  reflectance  of  the 
ceiling  cavity  at  the  upper  dashed  line  in  Figure  6.  This  converts  the  room 
interior  into  a  rectangular  parallepiped  of  height  h,  to  which  the  customary 
zonal  cavity  method  of  interior  lighting  design  may  be  applied.  The  room  cavity 
ratio  is  calculated  and  entered  in  location  9  on  Worksheet  2.  Then  the  room 
coefficient  of  utilization  is  obtained  from  Table  8  on  page  34,  a  table 
to  be  developed  in  the  next  portion  of  this  report.  The  result  is  entered  in 
location  10. 

Finally,  the  last  step  on  Worksheet  2  is  the  entry  of  the  expected  room 
surface  dirt  depreciation  (RSDD). 
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The  equations  in  steps  7  and  8  on  Worksheet  2  may  be  justified  as  follows: 
The  weighted  average  reflectance  of  the  ceiling  surface  is 


av 


pcAc+psAs 

VAs 


(19) 


where  A$  is  the  skylight  area  and  Ac  is  the  area  of  the  nonskylight  portion  of 


the  ceiling.  A$  and  Ac  are  given  by 


A$  =  2nab 

A  =  -2nab 
c  cos.|> 

Substituting  Equations  20  into  Equation  19  yields 


(20) 


p  =  p 

av  ^c 


2nab  cosa  ,  x 
- m - *  (pc-ps} 


(21) 


Equation  21  ignores  the  triangular  ends  of  the  ceiling  cavity,  which  are 
usually  relatively  small. 

With  oav  calculated,  we  can  proceed  to  obtain  pcc.  The  general  form  for 
the  effective  cavity  reflectance  of  a  non-horizontal  ceiling  is 

1 


cc 


£l  i-pav 
a„  ( — — ) 


(22) 


av 


where  A  is  the  area  of  the  ceiling  cavity  and  Aq  is  the  area  of  the  ceiling 
cavity  opening.  From  Equation  20,  Equation  22  is  rewritten  as 


1 


cc 


1- 


1  1-p 


(23) 


cos4>  (- 


av 


■) 


av 


A  comment  on  the  two  dirt  depreciation  factors  is  in  order.  SDD  represents 
dirt  accumulation  on  the  exterior  and  interior  of  the  skylights,  which  will 


reduce  the  transmittance  :  over  time.  The  authors  are  not  aware  of  any  data  on 
this  factor.  It  is  suggested  that,  in  the  absence  of  such  data,  a  value  of  .75  be 
assumed  and  that  the  skylights  be  cleaned  at  least  once  annually. 
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To  determine  RSDD,  the  method  in  IES  Handbook  VI,  Reference  Volume,  Chapter  9 
may  be  used,  with  the  skylight  interpreted  as  a  direct  luminaire.  Otherwise 
an  average  value  of  .90  may  be  assumed. 

We  turn  next  to  Worksheet  3,  which  requires  the  entry  of  sky  data.  The 
first  four  items  are  self-explanatory.  For  now,  the  sky  condition  is  either 
overcast  or  clear,  but,  in  the  future,  data  for  other  sky  conditions,  such  as 
partly  cloudy  may  become  available. 

Item  5  on  Worksheet  3  asks  that  the  illuminance  on  the  tilted  skylight  be 
entered  for  either  of  two  building  orientations.  If  the  walls  labeled  "i"  in 
Figure  4  face  north  and  south,  the  illuminance  on  the  north-facing  skylights  is 
entered  in  the  slot  and  that  on  the  south-facing  skylights  in  the  E^  slot. 

If  the  "i"  walls  face  east  and  west,  E^  is  the  illuminance  on  the  east-facing 
skylights  and  E ^  is  the  illuminance  on  those  skylights  facing  west.  The 
values  for  E1  and  E^  are  taken  from  the  appropriate  tables  in  Section  B  of  this 
report. 

All  information  for  the  calculation  of  average  maintained  interior 
illuminance  on  the  work  plane  has  now  been  assembled  and  we  turn  to  Worksheet  4. 

The  work  plane  area  and  skylight  area  per  side  are  entered.  Next  the  overall 
utilization  coefficient,  which  is  the  product  of  the  skylight  transmittance  and  the 
room  coefficient  of  utilization,  is  calculated  and  entered.  Following  this  the 
overall  light  loss  factor,  including  both  skylight  and  room  surface  dirt  de¬ 
preciation  is  calculated  and  entered.  Finally,  the  desired  illuminance  is 
obtained.  The  quantity  ( E-j  +E2 )  A$  is  the  luminous  flux  impinging  on  the  skylights. 
(E-j+E^A^  is  the  luminous  flux  impinging  on  the  work  plane  if  all  surfaces  are 
clean.  Multiplying  this  flux  by  gives  the  average  luminous  flux  over  time 
reaching  the  work  plane.  Division  by  the  work  plane  area  gives  the  average 
maintained  work  plane  illuminance. 
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A  table  of  room  coefficients  of  utilization  for  skylights  is  presented 
in  Figure  9-76  of  IES  Handbook  VI,  Reference  Volume.  These  are  given  for 
ceiling  reflectances  of  75%  and  50%  and  wall  reflectances  of  50%  and  30%,  as 
a  function  of  a  quantity  called  Room  Ratio  (RR).  These  data  may  be  made  a 
function  of  room  cavity  ratio  through  the  relation 

RCR  =  ^  (24) 

If  the  data  are  plotted  against  RCR  for  the  various  reflectances,  four  curves 
result  which  are  very  close  to  each  other  and  may  be  approximated  by  a  single 
straight  line  equation  as 

RCU  =  -.057RCR+.965  (25) 

Use  of  this  equation  eliminates  the  need  for  calculating  p3  ,  and  in  Steps 

flV  CC 

7  and  8  of  Worksheet  2. 

RCU's  from  Equation  25,  when  compared  with  those  in  Figure  9-76  in  the  IES 
Handbook,  differ  by  a  maximum  of  6%  in  the  range  1<RCR<8  and  by  slightly  higher 
percentages  outside  this  RCR  range. 
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FIGURE  6 


WORKSHEET  1 


Length 

Width 

Number  Per  Side 

Dirt  Depreciation  Factor 

Transmittance 


-  SKYLIGHTS 


Reflectance 


WORKSHEET  2  -  ROOM 


Directions  Roof  Faces 

Roof  Pitch 

Work  Plane  Length 

Work  Plane  Width 

Room  Cavity  Height 

Room  Surface  Reflectances 


N-S 
♦  ■ 
i  = 
w  = 
h  = 


pc  = 


pw  = 


Average  Ceiling  Reflectance 


pf  = 


2nab  cos*  ,  v 
- ZTo -  VP^-Pc)  = 


pav  ~  pc  wt  '"c  "s' 


Effective  Ceiling  Cavity  Reflectance 
1 


cc 


1  + 


1  1-p 


cos*  (- 


av 


av 


Room  Cavity  Ratio 


RCR  =  = 

wt 


Room  Coefficient  of  Utilization  RCU  = 
Dirt  Depreciation  Factor  RSDD  = 


,  E-W 
deg. 
ft. 
ft. 
ft. 
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WORKSHEET  3  -  SKY  ILLUMINANCE 


1. 

2. 

3. 

4. 

5. 


Oate 

Time 

Sky  Condition 
Lati tude 

Skylight  Illuminance 

E,  =  E  or  Ea 
I  n  e 

E,  =  Ec  or  E 
c  s  w 


E,  =  _  ft-c 

E2  =  _ ft-c 


WORKSHEET  4  -  INTERIOR  ILLUMINANCE 


1. 

Work  Plane  Area 

A  „  -  'Hi  - 
wp  - 

2. 

Skylight  Area  Per  Side 

A$  =  nab  =  _ 

3. 

Utilization  Coefficient 

Ku  =  RCUxt  =  _ 

4. 

Light  Loss  Factor 

Km  =  SDDxRSDD 

5. 

Average  Maintained  Work  Plane 

Illuminance 

_  (E,+E2)xKu«KmxAs 

wp  A 

r  wp 

ft-c 
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D.  ROOM  COEFFICIENTS  OF  UTILIZATION  FOR  SKYLIGHTS  IN  HORIZONTAL  ROOFS 


In  this  section  a  suggested  revision  of  Table  2  in  Navy  Tech  Data  Sheet 
80-09  will  be  presented.  New  values  for  the  room  coefficients  of  utilization 
appearing  in  this  table  will  be  calculated  using  the  procedure  in  IES  Handbook 
VI,  Reference  Volume,  pages  9-37  to  9-39.  This  procedure  may  be  adapted  to 
skylights  as  follows: 

1.  For  each  skylight,  define  nine  10-degree  flux  zones.  Let 
n  be  the  zone  index,  where  l<n<9,  with  zone  1  adjacent 

to  the  nadir  and  zone  9  adjacent  to  the  roof. 

2.  Compute  the  lumens  in  each  zone  from 


3. 


4. 


4>n  =  2irl(cose2-cos0^ ) 

where  and  e-|  are  the  angles  bounding  the  10-degree 
zone  and  I  is  the  mid-zone  intensity. 

Determine  the  total  (downward)  lumens  leaving  the  skylight. 

*  ■  X  ♦» 

Separate  the  total  lumens  into  two  parts,  those  traveling 
directly  to  the  walls  and  those  going  directly  to  the  work 
plane.  The  fraction  of  the  flux  which  reaches  the  work 
plane  directly  is 


D  =  -  7  k  * 
m  <p  rwrn 

where  k  is  called  a  zonal  multiplier  and  is  a  function  of 
mn 

cavity  ratio  (m),  zone  (n),  and  spacing  to  mounting  height 
ratio  (S/MH).  Zonal  multipliers  have  been  published  only  for 
S/MH  values  of  .4,  .7  and  1.0.  In  skylight  installations, 
the  S/MH  ratio  is  often  larger.  With  the  aid  of  unpublished 


(26) 


(27) 


(28) 
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data  assembled  by  R.E.  Levin  (GTE  Sylvania),  we  have  developed 
a  set  of  zonal  multipliers  for  an  S/MH  ratio  of  1.5.  These 
appear  in  Table  7  and  are  valid  for  room  length  to  width  ratios 
between  1  and  2  and  for  an  edge  ratio  of  1/2  (distance  from 
skylight  to  wall  is  one-half  distance  between  skylights). 


TABLE  7  ZONAL  MULTIPLIERS,  S/MH  =1.5 


Zone 


(n) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

2 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

3 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

4 

.98 

.96 

.94 

.92 

.90 

.88 

.86 

.84 

.82 

.80 

5 

.87 

.74 

.60 

.47 

.34 

.21 

.08 

0 

0 

0 

6 

.82 

.64 

.48 

.32 

.18 

.06 

0 

0 

0 

0 

7 

.75 

.52 

.33 

.18 

.09 

.01 

0 

0 

0 

0 

8 

.56 

.22 

.04 

0 

0 

0 

0 

0 

0 

0 

9 

.10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5. 

Obtain  the 

form 

factor  (f). 

Form 

factors  for 

various  room 

cavity  ratios,  based  on  L/W  =  1.6,  are  given  in  Figure  9-17 
in  the  IES  Handbook. 

Calculate  the  room  coefficient  of  utilization  from 

.  i-S'W'-V  . 

CU  "  m(l-Pl)(l-p3)C0 


p3c3(c1+c2) 

L  ih'-rs  V  J 


D_ 


n-p3)c0  J  i-p3 


(29) 
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-  j  -  a  ...  . 


where 


C,  = 


2.5p1(1-r)+mf(l-p1) 


0-p2)(H0 


i+P2f 


(30) 


C,  » 


(l-p3)(l+f) 


l+p3f 


Co  =  C1  +  C2  +  C3 


and  p-j ,  and  p^  are  the  reflectances  of  the  walls,  ceiling 
and  floor,  respectively. 

In  carrying  out  this  procedure  to  determine  room  coefficients  of 
utilization  for  skylights,  two  normalized  intensity  distributions  were  assumed 
namely. 


I  =  cose,  I  =  cos  9 


(31) 


The  first  represents  a  perfectly  diffusing  skylight,  whereas  the  second  is  a 
more  directed  distribution.  The  resulting  coefficients  of  utilization  appear 
in  Tables  8  and  9. 

An  examination  of  Tables  8  and  9  shows  that  the  coefficients  are  con- 

3  . 

sistently  greater  for  I  =  cos  0,  for  a  given  RCR,  pc  and  p  .  This  is  because 

-i 

of  the  greater  directivity  of  the  cos  0  distribution,  i.e.  more  flux 
reaching  the  work  plane  directly. 

It  should  not  be  assumed  that  either  of  the  intensity  distributions  in 
Equations  31  is  typical  of  a  particular  skylight.  Until  intensity 
distribution  curves  for  a  variety  of  skylights  under  a  variety  of  sun  angle 
and  atmospheric  conditions  become  available,  it  is  fruitless  to  speculate  on 
what  "typical"  distributions  might  be.  In  the  interim,  it  is  suggested  that 
the  cosy  (perfectly  diffusing)  distribution  be  used  as  the  norm  and  that 
coefficients  for  that  distribution  be  used  (Table  8).  Parenthetically,  it 
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should  be  noted  that  the  coefficients  in  Table  8  are  generally  larger  than  the 
corresponding  coefficients  obtained  from  Figure  9-76  in  the  IES  Handbook, 
which  have  been  used  in  skylight  calculations  for  a  considerable  number  of 
years. 
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pc 

RCR 

pw 

50% 

30% 

10%  | 

c* 

o 

00 

0 

1.19 

1.19 

1.19 

1 

1.11 

1.08 

1.06 

2 

1.03 

0.99 

0.96 

3 

0.92 

0.88 

4 

0.91 

0.86 

0.82 

5 

0.86 

0.80 

0.76 

6 

0.81 

0.75 

0.71 

7 

0.76 

0.71 

0.67 

8 

0.74 

0.68 

0.65 

9 

0.73 

0.68 

0.64 

10 

0.72 

0.67 

0.63 

i 

0 

1.11 

1.11 

i.n 

1 

1.05 

1.03 

1.01 

2 

0.98 

0.95 

0.93  1 

3 

0.93 

0.89 

0.76 

4 

0.88 

0.84 

0.80 

5 

0.83 

0.79 

0.75 

6 

0.79 

0.74 

0.70 

7 

0.74 

0.70 

0.66 

8 

0.72 

0.68 

0.64 

9 

0.71 

0.67 

0.64 

10 

0.71 

0.66 

0.63 

20% 

0 

1.04 

1.04 

1.04  | 

1 

1 

1 

0.99 

0.98 

0.97  I 

| 

2 

0.94 

0.92 

| 

3 

0.89 

0.87 

0.84 

i 

4 

0.85 

0.82 

0.79 

5 

0.81 

0.77 

0.74 

6 

0.73 

7 

0.73 

0.69 

0.66 

8 

0.71 

0.67 

0.64  1 

9 

■5E— 

0.66 

10 

KM 

0.66 

0.63 

TABLE  9  COEFFICIENTS  OF  UTILIZATION 
pf  =  20%,  S/MH  =  1.5,  I  =  cos 30 
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£.  APPENDIX 


A.  Derivation  of  overcast  sky  equations  for  illuminance, 

B.  Derivation  of 

C.  Derivation  of  (Ejsj).  total  solar  illuminance. 

D.  Derivation  of  clear  sky  equations  for  illuminance. 
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A.  Derivation  of  Overcast  Sky  Equations  for  Illuminance  Due  to  Sky  luminance 


1.  The  horizontal  illuminance  for  a  sloped  surface  under  overcast  sky 
conditions  is  given  by  Equation  Al. 

EH0C  =  EHFLAT  *  EHSTRIP  (A1) 

where  EHFLAT  =  horizontal  illuminance  due  to  overcast  sky  luminance 

of  entire  sky  dome. 

EHSTRIP  =  hor i zontal  illuminance  due  to  overcast  sky  luminance  of 
a  strip  of  the  sky  dome  (See  Figure  Al ) 


Figure  Al 


2.  The  vertical  illuminance  is  due  to  the  uniform  luminance  distribution  with 
respect  to  <t> .  Each  e-ring  has  constant  luminance  all  the  way  around.  It  is 
seen  that  the  only  vertical  component  is  due  to  the  strip's  non-contribution 
(if  the  sky  dome  is  divided  into  half  with  the  left  of  the  surface  as  positive 
and  right  of  the  surface  negative)  and  is  positive. 

,r  _  LcoscdA 

atv  "  2 

r 

where  cose  =  sinesin<j> 

dEv  X 1  +2cose ) s in^si n  eded4> 

r*-<P'  fTt/2  L  2 

E  =  sin<j>  -4  (l+2cose)sin  eded4> 

'  <, p'  7r/2-a 

which  is  approximated  by 
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C.  Derivation  of  Total  Solar  Illuminance 

By  approximation  of  the  sun  as  a  point  source  and  using  the  inverse  square 
law 

c  -  lCOS> 


where 


But 


y  =  angular  distance  between  sun  and  surface  normal  and  is  given  by 


-TSI 


•PSI 


'TSI 


-HSI 


'TSI 


y  =  cos"^ (sinSAcosa+cosSAsinacos(B-SAZ)) 
(Reference  5) 

=  -i  (sinSAcosut+cosSAsinacos(6-SAZ) 
d^ 

J_ 

"  d2 

=  EpSIsinSAcosa+Ep^jCOsSAsinacos(B-SAZ) 

=  EpSjSinSA  and  EVSJ  =  EpSjCOsSA 

=  EHSIcosa+EVSIsinaCOs(2-SAZ) 


D.  Derivation  of  Clear  Sky  Equations  for  Illuminance  Due  to  Sky  Luminance 
1.  As  in  the  overcast  sky  case  the  horizontal  illuminance  can  be  calculated  by 

ehc  =  ehcflat'  ehcstrip 

where  ^HCFLAT  =  honzontal  illuminance  due  to  clear  sky  luminance  of 

the  entire  sky  dome 

EHCSTRIP  =  hor1zonta^  illuminance  due  to  clear  sky  luminance  of 
a  strip  of  the  sky  dome  (see  Figure  A1 ) 

The  portion  of  horizontal  illuminance  for  each  sky  element  is 

L  cosodA 

AC  C 
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where  dA  =  (rsine)d4>rde 

L  =  Lz(1-e~‘32/coS^avg)(.91-»-10e'3^.45cos2Yl 
c  .279(.91+10e'3(7r/2  SA)  +  .45sin2SA 
The  clear  sky  luminance  (Lc)  is  an  approximation  where 


=  ARCcosfsinSAcose,  -cosSAsine.^sinU  -SAZ)] 
91  +  02 


'avg  =  — — 

pavg  =  — 2 —  (see  figure  4) 


The  horizontal  illuminance  on  a  flat  surface  is 


-HCFLAT 


Lcsinecoseded$ 


By  using  an  approximation  of  Lc  and  summing  by  taking  small  increments,  in 

first  the  e  then  the  $  direction,  the  above  integral  becomes 

2  2 

-jt/2  sin  e,-sin  e-,  2 n 

EHCFLAT  =  l  2  ^  l  L(*2"h)j’ 

6_0,einc  *'°’»inc 

When  the  surface  is  not  flat  the  surface's  orientation  becomes  a  factor 
because  the  contribution  of  that  part  of  the  sky  that  is  "beneath"  the  surface's 
plane  projection  to  the  sky  dome  is  lost.  The  clear  sky  luminance  is  not 
uniform  with  respect  to  the  <f>  direction  so  the  contribution  to  horizontal 
illuminance  will  vary  with  orientation  of  the  surface. 

The  horizontal  illuminance  due  to  the  strip  beneath  the  surface's  plane 
projection  that  intersects  the  sky  dome  (see  Figure  5)  is 


"HCSTRIP 


*-<$'+8  w/2 

l  l  L  sinecosededij) 

<j>  +  8  TT 

2  a 


where 


■f'  =  intersection  of  sky  dome  and  surface  plane  projection  for  each 
"ring"  of  sky  elements  (see  Figure  4) 

8  =  angle  of  orientation  with  respect  to  south  and  counter-clockwise 

direction  being  positive  (see  Figure  2) 
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The  approximate  horizontal  strip  illuminance  is 

2  2 

tt/2  sin  02-sin 

ehcstrip  =  ^  ^  ?  ) 


v-ip'+B 

I  L(4>2-Ih  )3 


°~2  'a,einc 


2.  The  vertical  illuminance  of  the  clear  sky  is 


EVC  =  EVCP0S  "  ^EVCNEG  ■  EVCSTRIp)- 

where  ^VCPOS  =  P0Sltlve  vertical  illuminance  component  (half  of  sky  dome 

without  strip). 

EVCNEG  =  ne9ative  vertical  illuminance  component  (half  of  sky 


dome  with  strip). 


EVCSTRIP  =  vertica^  illuminance  due  to  strip. 

The  portion  of  vertical  illuminance  for  each  sky  element  is 


dEy  = 


LccoscdA 


where  cos c  3  sinesin^, 

dA  =  (rsine)d<j>rde 


Lc  =  sky  luminance  of  clear  sky  (same  as  in  Appendix  D1 ) . 


The  positive  vertical  component  is  then 


2tt+S  2t 


:vcpos 


.  2 

Lcsin4>sm  eded4> 


r+B  0 


and  is  approximated  by 


tt/2  e9-0i  sin20,-sin202  2tt+b 

EVCP0S  =  £  (  2  +  4  ^  £  L  |  (cos<j>i -cosi^2 )  I  • 

e=0,9inc  *=1T+6’*inc 


The  negative  vertical  component  is 


■VCNEG 


n  r+b  tt  it 
)  ft  nJ 


Lcs in<j>s in  eded4> 


and  is  approximated  by 
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VCNEG  *  '  2 

0  °’6inc 

The  strip  vertical  component  is 


tt/2  O/j-9-i  sin20.  -sin2e0  n+6 

I  (4-!-  *  1  2 


)  l  L 1  (cos<f>1-cos<i>2) 


*=6’*inc 


tt/2  2 

^VCSTRIP  =  J  |  Lcsin<i>sin  6 d 6 d<^ 
ir/2-a 

and  is  approximated  by 


-VCSTRIP 


tt/2  e9-e,  sin29,-sin20o  tt-a'+B 

l  (-V-  + - T - l 


L  |  ( cos4>-| -cos<j>2 ) 


9  =  2  ”a,0inc 


inc 


3.  The  total  illuminance  due  to  sky  luminance  of  the  completely  clear  sky  is 


Ec  =  E^c  cosa  +  E^  sina 
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NAVSE.C  CiRl  ACT  CO  (PWD).  Adak.  AK:  CC)  (Code  40U).  Ed/ell.  Scotland:  CO.  Cialeta  Ishnd.  Panama 
Canal.  CO  (Code  N60).  Homestead.  FI.:  CC)  (Code  30).  Puerto  Rico:  CO  (Energy  Conscrv).  Sonoma.  CA; 
CO  (Energy  Conserv.)  Winter  llarhor.  ME:  CikIc  40.  C  he- a|uake.  VA:  PWO.  Torri  Sta.  Okinawa 
NAVSECORCCOM  Energy  Conserv  .  Washington  DC 
NAVSECSTA  C  ode  540.  Washington  DC 

NAVSHIPYD  C  ode  402.4.  Philadelphia  PA:  CO  ( C  ode  405):  PW1)  (Code  400  03).  Charleston  SC:  Commander 
(Code  406).  Portsmouth.  Nil 
NA\  REC.MFDCEN  CO  (Code  133).  Long  Beach.  CA 
NAVSHIPYD  Puget  Sound.  CMDR  (Code  402  3).  Bremerton.  WA 

NAVS1A  CO  (PWD).  Kellavik.  Iceland:  CO  (PWD).  Rota.  Spam:  CO  (Code  18410).  Mayp.  .  FL;  CO  (Code 
413).  Grmo.  Cuba;  (Code  50A)  Rodman.  Panama  C  anal:  CO  (Code  ODE).  San  Diego.  CA.  CO  (Energy 
Conserv)  CC)  (Code  52).  Brooklyn  NY 

NAVSl'BASE  CO  (Code  803).  Groton.  CT:  CO  iCodc  223)  Bangor.  Bremerton.  WA 
NAVSl'PPACT  CO  (Code  N52).  New  Orleans.  LA:  CO  (Code  413).  Seattle.  WA:  CO  (Energy  Conserv). 
Naples.  Icily 

NAVSl'PPBASE  CO  (Energy  Conserv)  Kings  Bay.  CiA 

NAYSCPPFAC  Code  02.  Thurmonl  Ml):  CO  (Energy  Conserv)  Diego  Garcia  I 

VAVSCPPO  CO  (APWO)  la  Maddalena.  Italy 

NAVSC  REW PNC  IN  Dahlgren  lab.  WWH2  Dahlgren  VA 

NAVI'E.l  ('OMNICOM  C  ode  05.  Washington  DC 

NAVI  SI  AWAR!  NGSI  A  CO  1C  ode  070  2).  Kevpoil.  W  \ 

NAVWPNC  EN  Commander  (Code  2o>5|  (Irma  lake.  C  \ 

NANWI’NSEA  CC)  (Energy  Conserv).  Colls  Neek.  N.L  CO  (Code  0*)22ll.  C  oneorvl.  CA;  CO  (Energy  Conserv) 
Yorktown.  VA:  Code  (Nil.  Seal  lieaeh  CA 
NAVWPNSl  PPC  EN  CO  (Code  0021.1.  Crane.  IN 

NC  BC  CO  (Code  80).  Pori  Hucncme.  C  A.  CO  (Energy  C  onserv  I.  Davisville.  Rl 
NOAA  Library  Rockville.  Ml) 

NRl.  PWO  C  ode  2530.1.  Washington.  DC 

NSC  CO  (C  ode  70A).  Puget  Sound.  W  A.  CO  (C  ode  46.A)  San  Diego.  CA 
NSD  CO  (Cixle  50E) 

Nrt  CO  (C  ode  NAC50FI  Orlando.  1 1 

NL'SC  CC)  (Code  5204).  Newport.  Rl 

ONR  CO  (Code  701)  Pasadena.  CA 

PAC  MISRANFAC  CO  (Code  7031).  Kekuha.  III 

PMTC  Commander  (Code  6200-3).  Point  MuguL.  C  A 

PWC  CO  (Code  610).  Pensacola.  FI  ;  C  ode  KKIA.  Cireal  I  .ikes.  II 

PWC-NAS  NAS  Pensacola.  EE 

PWC  CO  (C  ode  613).  San  Diego.  CA:  C  ode  600A  Norfolk.  VA:  CO  1C  ode  10031.  Oakland.  CA:  CO  (Code 
153).  Guam:  CO  (Code  1001:3).  Oakland.  CA;  CO  (Code  )0|.  IVarl  Harbor.  Ill:  CO  (Code  IIIOE).  San 
Diego.  CA:  Code  116,  Seattle.  WA:  CO  (Code  601 1.  Suhie  Bav 
SKC  CO  (C  ode  763).  Meehanicsburg.  PA 
SCPSHIP  ADMINO.  San  Francisco.  CA:  CikIc  ‘ml 
I'SNA  C  ode  1 7(1.  Annapolis.  MD 
ARIZONA  Kroelinger  Fempe.  A/ 

AL'Bl'RN  F  NI\  Bldg  Set  Dept.  I.cchncr.  Auburn  \l 


C  i  \  MSI  IN  I'NIV  lol  Arch  .  I  g.m.  (  IvuNm.  Si 
IKANKMN  INiS Illl'll  M  KuIumv  PhiUlcIphia  PA 

OlORUIA  INS  Mil  II  Of  If rilSOl  (HiV  Col  Arch  Benton  Atlanta.  (i\ 
HARVARD  I'NIV  Dept  of  Aichitcciorc.  I>i  kirn.  Cambridge.  MA 
U )  V\  A  SIAM  UNIVERSITY  Dept  Arch.  McKroxxn.  Antes.  IA 
LAWRENC  E  BERK  L  AB  Window  A  l  ighting  Prog  Berkelcx .  c  \ 

I  .OS  ALAMOS  SC  I  l.AB  Solar  i  nergx  < ip.  I  os  Alamos.  NM 
MIL  Cambridge  MA  (Rm  10*501).  lech  Reports.  Engr  Uh  l 
PENNSYLVANIA  STATE  I'NIVFRSMY  Arch  Lngr  Dept.  College  Station.  PA 
I’NIV  COLORADO  Arch  Engr  Dept 

CNIVFRSITY  OF  ARIZONA  S  I)  Lakeman.  (C  ollege  ol  Architecture)  I  ucson  A/ 
VNIVERSITY  OF  FLORIDA  Dept  Arch  .  Morgan,  Ciumcsxillc.  LI. 

UNIVERSITY  OF  NEW  HAMPSHIRE  f  lee.  Lngr  Depot.  Di  Mimloch.  Durham.  N  II 

CNIVLRSITY  OF  WASHINGTON  Insti  lor  I'nvir  Studies 

PGA  I  lihrarx .  San  Francisco.  C  A 

SANDIA  LABORATORIES  Library  Dix.  Livermore  C  A 
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